From Spin Ladders to the 2-d 0(3) Model at Non-Zero Density 

S. Chandrasekharan" , B. Scarlet^ and U.-J. Wiese''''^ 

'^Department of Physics, Duke University, Box 90305, Durham, North Carolina 27708 
^ Center for Theoretical Physics, Laboratory for Nuclear Science and Department of Physics, 

Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 
'^Institute for Theoretical Physics, Bern University, Sidlerstrasse 5, 3012 Bern, Switzerland 



Abstract 

The numerical simulation of various field theories at non-zero chemical potential suffers from severe complex 
action problems. In particular, QCD at non-zero quark density can presently not be simulated for that reason. 
A similar complex action problem arises in the 2-d 0(3) model — a toy model for QCD. Here we construct the 
2-d 0(3) model at non-zero density via dimensional reduction of an antiferromagnetic quantum spin ladder in a 
magnetic field. The complex action problem of the 2-d 0(3) model manifests itself as a sign problem of the ladder 
system. This sign problem is solved completely with a meron-cluster algorithm. 
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1. Introduction 

Numerical simulations of numerous quantum 
systems suffer from notorious sign and complex 
action problems. For such systems the Boltzmann 
factor of a configuration in the path integral is in 
general complex and can hence not be interpreted 
as a probability. When the complex phase of the 
Boltzmann factor is incorporated in measured ob- 
servables, the fluctuations in the phase give rise 
to dramatic cancellations. In particular, for large 
systems at low temperatures this leads to relative 
statistical errors that are exponentially large in 
both the volume and the inverse temperature. As 
a consequence, it is impossible in practice to study 
such systems with standard importance sampling 
Monte Carlo methods. 

Recently, some severe sign problems have 



been solved with meron-cluster algorithms [1~5]. 
Meron-clusters are used to identify canceling pairs 
of configurations with the same weight but op- 
posite signs. Configurations with merons exactly 
cancel in the path integral. The Monte Carlo sim- 
ulation can thus be restricted to the zero-meron 
sector with positive weights for which standard 
importance sampling works efficiently. 

It is natural to ask if the meron-concept can be 
applied to the complex action problem in dense 
QCD. This is indeed the case in the limit of in- 
finitely heavy quarks in the Potts model approx- 
imation to QCD [6]. For dynamical light quarks, 
on the other hand, it is not obvious if the meron 
concept can be appHed to QCD. This seems most 
likely in the D-theory formulation of field theory 
in which 4-d QCD arises via dimensional reduction 
from a (4 + l)-d quantimr link model [7,8]. Phys- 
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ical gluons emerge as collective excitations of dis- 
crete variables [7] so-called quantum links — 
which are gauge covariant generalizations of quan- 
tum spins. In this paper we study the D-theory 
formulation of the 2-d 0(3) model — a toy model 
for QCD — which arises via dimensional reduction 
from an antiferromagnetic spin ladder. In this case, 
the discrete variables are ordinary quantum spins. 
At the level of the quantum spin system, the chem- 
ical potential of the 0(3) model manifests itself as 
an external magnetic field. 



2. Spin Ladders in a Magnetic Field and 
the 2-d 0(3) Model at Non-Zero Density 

Wc consider a ladder system of quantum spins 
1/2 on a square lattice of size LxL' {L' even, L ^ 
L') with periodic boundary conditions in both di- 
rections. The spins located at the sites x are de- 
scribed by operators S^; with the usual commuta- 
tion relations 

[^l^^i] = 'i'^xy^ijkS^. (1) 

The antiferromagnetic Hamilton operator (J > 0) 
if = J^S.-S,+j-B.^S„ (2) 

couples the spins at the lattice sites x and x + i, 
where i is a unit-vector in the i-direction. 

Chakravarty, Halperin and Nelson used a (2-1-1)- 
d effective field theory to describe the low-energy 
dynamics of spatially 2-d quantum antiferromag- 
nets [9]. Chakravarty has applied this theory to 
quantum spin ladders with a large even number 
of coupled spin 1/2 chains [10]. These systems are 
described by the effective action 

P L L' 

^[e] = J dt J dx J dy ^[dxG ■ dxG + dye ■ dyS 





+ \dte-dte\, 



(3) 



where ps is the spin stiffness and c is the spinwave 
velocity. 

When the spin ladder is placed in a uniform 
external magnetic field B, the field couples to a 



conserved quantity — the total spin. Hence, on 
the level of the effective theory, the magnetic field 
plays the role of a chemical potential, i.e. it appears 
as the time-component of an imaginary constant 
vector potential. As a consequence, the ordinary 
derivative dte is replaced by the covariant deriva- 
tive dte -|- iB X e and the action takes the form 

L L' 

S[e]= J dt J dx J dy ^[dxG ■ dxG + dyB ■ dyB 



+ ^{dte + iBxe)-{dte + iBxe)]. (4) 

For a sufficiently large number of coupled chains 
{L' ^ c/ps) the system undergoes dimensional re- 
duction to the 2-d 0(3) model with the action 

'S'N = j dt j dx ^^^[d^e ■ dxB 



+ ^{dte -I- iB X e) • {dte + iBx e)] 

/3c 



: j d{ct) j dx -^[dxG ■ dx 





+ {dcte + in X e) ■ {dcte + ip X e)]. (5) 

The effective coupling constant is given by l/g^ = 
PsL' / c and the magnetic field appears as a chemical 
potential of magnitude p = B/c. 



3. Path Integral for Quantum Magnets 

To derive a path integral representation of the 
partition function wc decompose the Hamilton op- 
erator of eq.(2) into five terms 



H ^ H1+H2 + ... + H5. 

The various terms take the form 



Hi — 2, ^x,i^ -Hj+2 — 



^ = (^1,^2) 
Xieven 



hx,i^ 



■ = (^1 >^2) 

Xiodd 



with hx,i = JSx ■ S^_^i and 



=(xi,X2) 



(6) 



(7) 



(8) 
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with hx = —BS].. The individual contributions 
to a given Hi commute with each other, but two 
different Hi do not commute. Using the Trotter- 
Suzuki formula we express the partition function 
as 



Z = Tr[exp(-/3il)] 

iVl — >OG 

xexp(-eif2)-exp(-eif5)]^ 



limTr[exp(-eifi) 

M— >oo 



(9) 



We have introduced 5M Euclidean time slices with 
e = (}/M as the lattice spacing in the Euclidean 
time direction. We insert complete sets of eigen- 
statcs I t) and | J,) with eigenvalues = ±1/2 
between the factors exp{—eHi). 

The partition function is now expressed as a path 
integral 



Z = 5^Sign[s]exp(-S[s]), 



(10) 



over configurations of spins s{x,t) =T,i on a 
(2 + l)-dimensional space-time lattice of points 
The Boltzmann factor cxp(— S'[s]) is a 
product of space-time plaquette contributions 
exp{-5[s(a;, t), s{y, t), s{x,t+l), s{y,t + l)]} with 

exp(-5[t, T, T, T]) = exp(-5U, i, i, i]) = 

= exp(-eJ/2), 
exp(-5[T, i, T, i]) = exp(-5[i, T, i, T]) = 

= cosh(eJ/2), 
exp(-5[T, i, i, T]) = exp(-5[i, T, T, i]) = 

= sinh(eJ/2), (11) 

as well as the time-like bond contributions 

exp{— S'[s(x, t), s{x, t + 1)]} with 

cxp(-5[T, T]) = exp{-S[i, i]) = cosh(ei?/2), 
exp(-5[T, I]) = exp{-S[l, T]) = sinh(eS/2). (12) 

The sign of a configuration, Sign[s], also is a 
product of space-time plaquette contributions 
Sign[s(a;, t), s{y, t), s{x, t + 1), s{y, t + 1)] with 

Sign[T,T,T,T]) = Sign[i,i,i,i]) = l, 

Sign[t,i,T,i]) = Sign[i,t,i,T]) = l, 

Sign[t, i, i, T]) = Sign[i, T, T, I]) = -1. (13) 

Figure 1 shows two spin configurations in (1 -I- 1) 
dimensions. The first configuration is completely 



antiferromagnetically ordered and has Sign[s] = 1. 
The second configuration contains one interaction 
plaquette with configuration [J,, T) Ti i] which con- 
tributes Sign[J,, t, t, I] = —1, such that the whole 
configuration has Sign[s] = — 1. 

The central observable of our study is the uni- 
form magnetization M = J2x 
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Fig. 1. Two spin configurations in (1 + 1) dimensions. The 
shaded plaquettes and tviae-like bonds carry the interac- 
tion. Filled dots represent spin up and open circles repre- 
sent spin down. The first configuration has Sign{s\ = 1 and 
the second configuration has Sign[s\ = — 1. The fat black 
line represents a meron- cluster. The other clusters are not 
shown. Flipping the meron- cluster changes one configura- 
tion into the other and changes Sign[s\ . 



4. Meron-Cluster Algorithm 

The meron-cluster algorithm is based on a clus- 
ter algorithm for a modified model without the 
sign factor. Quantum spin systems without a sign 
problem can be simulated very efficiently with the 
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loop-cluster algorithm [11-14]. The idea behind 

the algorithm is to decompose a configuration into 
clusters which can be flipped independently. Each 
lattice site belongs to exactly one cluster. When 
the cluster is flipped, the spins at all the sites on 
the cluster are changed from up to down and vice 
versa. The decomposition of the lattice into clus- 
ters results from connecting neighboring sites on 
each space-time interaction plaqiicttc or time-like 
interaction bond according to probabilistic cluster 
rules. A set of connected sites defines a cluster. In 
this case the clusters are open or closed strings. 
The cluster rules are constructed so as to obey de- 
tailed balance. To show this property we first write 
the space-time plaquctte Boltzmami factors as 

exp(-S'[s(a:;, t),s{y, t), s{x, t + 1), s{y, t + 1)]) = 

■^^s{x,t),s(x,t+l)'^s{y,t),s{y,t+l) + 

B5s(x,t),-s{yJ)(^s{x,t+l),-s{yJ + l)- (14) 

The 5-functions specify which sites are connected 
and thus belong to the same cluster. The coefH- 
cients A and B determine the relative probabilities 
for different cluster break-ups of an interaction pla- 
quette. For example, A determines the probabil- 
ity with which sites are connected with their time- 
like neighbors and B determines the probability 
for connections with space-like neighbors. Insert- 
ing the expressions from eq.(ll) one finds 

exp(-5[T, T, T, T]) = exp(-5[i, i, j, i]) = 

= exp(-eJ/2) = A, 
exp(-S[T, i, T, i]) = cxp(-5[i, T, i, T]) = 

= cosh(gJ/2) = A + B, 
exp(-5[T, i, I, T]) = exp(-5[i, T, T, i]) = 

= sinh(e J/2) = S. (15) 

Similarly, the time-like bond Boltzmann factors are 
expressed as 

exp(— S'[s(x, t), s{x, t + 1)]) = 

C5s(^x,t),s{x,t+i) + D. (16) 

The probability to connect spins with their time- 
like neighbors is C /{C + D). The spins remain dis- 
connected with probabihty D/{C -\- D). Inserting 
the expressions from eq.(12) one obtains 



exp(-5[T,T]) = exp(-5[i,i]) = 

= cosh(eB/2) =C + £), 

exp(-5[T,i]) = exp(-SU,T]) = 

= sinh(eB/2) = D. (17) 



The cluster rules are illustrated in table 1. 
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Table 1 



Cluster break-ups of various plaquette and time-like bond 
configurations together with their relative probabilities 
A,B,C',D. Filled dots represent spin up, open circles rep- 
resent spin down, and the fat lines are the cluster connec- 
tions. 

Eqs.(14,16) can be viewed as a representation 
of the original model as a random cluster model. 
The cluster algorithm operates in two steps. First, 
a cluster break-up is chosen for each space-time in- 
teraction plaquette or time-like interaction bond 
according to the above probabilities. This effec- 
tively replaces the original Boltzmann weight of 
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a configuration with a set of constraints repre- 
sented by the (^-functions associated with the cho- 
sen break-ups. The constraints imply that the spins 
in one cluster can only be fiipped together. Sec- 
ond, every cluster is flipped with probability 1/2. 
When a cluster is flipped, the spins on all sites that 
belong to the cluster are flipped from up to down 
and vice versa. Eqs.(15,17) ensure that the cluster 
algorithm obeys detailed balance. 

The above cluster rules were first used in a sim- 
ulation of the Hciscnbcrg antifcrromagnct [12] in 
the absence of a magnetic field. In that case there 
is no sign problem. Then the corresponding loop- 
cluster algorithm is extremely efficient and has al- 
most no detectable autocorrelations. When a mag- 
netic field is switched on the situation changes. 
When the magnetic field points in the direction of 
the spin quantization axis (the 3-direction in our 
case) there is no sign problem. However, the mag- 
netic field then explicitly breaks the Z(2) fiip sym- 
metry on which the cluster algorithm is based, and 
the clusters can no longer be flipped with probabil- 
ity 1/2. Instead the flip probability is determined 
by the value of the magnetic field and by the mag- 
netization of the cluster. When the field is strong, 
flips of magnetized clusters are rarely possible and 
the algorithm becomes inefficient. To avoid this, 
we have chosen the magnetic field to point in the 
1-direction, i.e. perpendicular to the spin quantiza- 
tion axis. In that case, the cluster fiip symmetry is 
not affected by the magnetic field, and the clusters 
can still be flipped with probability 1/2. However, 
one now faces a sign problem and the cluster algo- 
rithm becomes extremely inefficient again. Fortu- 
nately, using the meron concept the sign problem 
can be eliminated completely and the efficiency of 
the original cluster algorithm can be maintained 
even in the presence of a magnetic field. 



5. Meron-CIusters and the Sign Problem 

Let us consider the effect of a cluster fiip on the 
sign. The flip of a meron-cluster changes Sign[s], 
while the flip of a non-meron-cluster leaves Sign[s] 
unchanged. An example of a meron-cluster is 
shown in figure 1. When the meron cluster is 



flipped, the first configuration with Sign[s] = 1 
turns into the second configuration with Sign[s] = 
— 1. This property of the cluster is independent 
of the orientation of any other cluster. Since fiip- 
ping all spins leaves Sign[s] unchanged, the total 
number of meron-clusters is always even. 

The meron concept allows us to gain an expo- 
nential factor in statistics. Since all clusters can 
be fiipped independently with probability 1/2, one 
can construct an improved estimator for (Sign) by 
averaging analytically over the 2^"^ configurations 
obtained by fiipping the Nq clusters in a configu- 
ration in all possible ways. For configurations that 
contain merons, the average Sign[s] is zero because 
fiipping a single meron-cluster leads to a cancel- 
lation of contributions ±1. Hence only the con- 
figurations without merons contribute to (Sign). 
The probability for having a configuration without 
merons is equal to (Sign) and is exponentially sup- 
pressed with the space-time volume. The vast ma- 
jority of configurations contains merons and con- 
tributes an exact to (Sign) instead of a statis- 
tical average of contributions ±1. In this way the 
improved estimator leads to an exponential gain 
in statistics. One can show that the contributions 
from the zero-meron sector are always positive. 

One can also find a simple expression for the 
improved estimator for the magnetization (M^) in 
terms of a winding number Wi of closed loops which 
result from joining open string clusters. One can 
define Wi for each loop to be its temporal wind- 
ing number. If a particular loop is not composed 
of open string clusters then Wi =0. With this def- 
inition of Wi it is easy to show that 

where iV is the number of meron-clusters. 

Since the magnetization gets non- vanishing con- 
tributions only from the zero-meron sector, it is un- 
necessary to generate any configuration that con- 
tains meron-clusters. This observation is the key 
to the solution of the sign problem. In fact, one 
can gain an exponential factor in statistics by re- 
stricting the simulation to the zero-meron sector, 
which represents an exponentially small fraction of 
the whole configuration space. We visit all plaque- 
tte and time-like bond interactions one after the 
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other and choose new pair connections between the 
sites according to the above chister rules. A newly 
proposed pair connection that takes us out of the 
zero-meron sector is rejected. After visiting all pla- 
quette and time-like bond interactions, each clus- 
ter is flipped with probability 1 /2 which completes 
one update sweep. In practice, it is advantageous 
to occasionally generate configurations containing 
merons even though they do not contribute to our 
observable, because this reduces the autocorrela- 
tion times. 

Figure 2 shows the magnetization density of an- 
tiferromagnetic quantum spin ladders with L' — i 
compared to analytic results (valid for large L') ob- 
tained with the Bethe ansatz [15] using the mass- 
gap m = 0.141(2)/a [16] and the spinwave velocity 
c = 1.657(2) Ja [17]. The agreement is remarkable 
and involves no adjustable free parameters. 




Fig. 2. Magnetization density {M^)/L of quantum spin lad- 
ders consisting of L' = 4 coupled chains as a function of 
the magnetic field B. The numerical data are for two sys- 
tems: one of size L = 20 at inverse temperature pj = 15 
(circles) and the other for L = 40 at fij = 24 (diamonds) . 
The solid curve is the infinite volume, zero temperature an- 
alytic result, uihile the two dashed curves are finite volume, 
non-zero temperature analytic results for the two simulated 
systems in the intermediate B region. The dotted curve 
represents saturation of the magnetization per spin at 1/2. 



6. Conclusions 

Using D-theory, the 2-d 0(3) model at non-zero 
chemical potential has been obtained from dimen- 



sional reduction of a (2 + l)-d quantum spin ladder 
in a magnetic field. The resulting sign problem has 
been solved completely with a meron-cluster algo- 
rithm. This is the first time that this toy model 
for QCD has been simulated efficiently at non-zero 
chemical potential. The next challenge is to address 
the complex action problem of dense QCD. In D- 
theory quarks arise as domain wall fermions and 
gluons emerge as collective excitations of quantum 
links. Hence, one needs to generalize the meron- 
concept to quantum link models as well as to do- 
main wall fermions. 

References 

[1] W. Bietenholz, A. Pochinsky and U.-J. Wiese, Phys. 
Rev. Lett. 75 (1995) 4524. 

[2] S. Chandrasekharan and U.-J. Wiese, Phys. Rev. Lett. 
83 (1999) 3116. 

[3] S. Chandrasekharan, J. Cox, K. Holland and U.-J. 
Wiese, Nucl. Phys. B576 (2000) 481. 

[4] J. Cox, C. Gattringer, K. Holland, B. Scarlet and U.-J. 
Wiese, Nucl. Phys. (Proc. Suppl.) 83 (2000) 777. 

[5] ^ Chandrasek haran and J. Osborn, 



:ond-mat/0109424 



[6] M. Alford, S. Chandrasekharan, J. Cox and U.-J. 
Wiese, Nucl. Phys. B602 (2001) 61. 

[7] S. Chandrasekharan and U.-J. Wiese, Nucl. Phys. 
B492 (1997) 455. 

[8] R. Brower, S. Chandrasekharan and U.-J. Wiese, Phys. 
Rev. D60 (1999) 094502. 

[9] S. Chakravarty, B. L Halperin and D. R. Nelson, Phys. 
Rev. Lett. 60 (1988) 1057; Phys. Rev. B39 (1989) 2344. 

[10] S. Chakravarty Phys. Rev. Lett. 77 (1996) 4446. 

[11] H. G. Evertz, G. Lana and M. Marcu, Phys. Rev. Lett. 
70 (1993) 875. 

[12] U.-J. Wiese and H.-P. Ying, Z. Phys. B93 (1994) 147. 

[13] H. G. Evertz, The loop algorithm, in Numerical 
Methods for Lattice Quantum Many-Body Problems, 
ed. D. J. Scalapino, Addison- Wesley Longman, 
Frontiers in Physics. 

[14] B. B. Beard and U.-J. Wiese, Phys. Rev. Lett. 77 
(1996) 5130. 

[15] P. B. Wiegmann, Phys. Lett. B152 (1985) 209. 

[16] O. F. Syljuasen, S. Chakravarty and M. Greven, Phys. 
Rev. Lett. 78 (1997) 4115. 

[17] B. B. Beard, R. J. Birgeneau, M. Greven and U.-J. 
Wiese, Phys. Rev. Lett. 80 (1998) 1742. 



6 



